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Abstract 

We study the structure of the normal matrix model (NMM). We show that all 
correlation functions of the model with axially symmetric potentials can be expressed in 
terms of holomorphic functions of one variable. This observation is used to demonstrate 
the exact solvability of the model. The two-point correlation function is calculated in 
the scaling limit by solving the BBGKYQ chain of equations. The answer is shown 
to be universal (i.e. potential independent up to a change of the scale). We then 
develop a two-dimensional free fermion formalism and construct a family of completely 
integrable hierarchies (which we call the extended--ftTP(iV) hierarchies) of non-linear 
differential equations. The well-known KP hierarchy is a lower-dimensional reduction 
of this family. The extended- KP{\) hierarchy contains the (2+l)-dimensional Burgers 
equations. The partition function of the (NxN) NMM is the r function of the extended- 
KP{N) hierarchy invariant with respect to a subalgebra of an algebra of all infinitesimal 
diffeomorphisms of the plane. 



1 Born, Bogoliubov, Green, Kirkwood, and Yvone, J32|. 
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1 Introduction 



Normal Matrix Model (NMM) was introduced in || and its connection to the Quantum 
Hall Effect was indicated. In || this connection was given a precise form by observing 
that the partition function of NMM coincides with zero-temperature partition function of 
two-dimensional electrons in the strong magnetic field. 

In the present paper we study the mathematical structure of the NMM. The distinct 
feature of the NMM, in contrast to other matrix models (Hermitian, Unitary, etc.), is its 
relation to two-plus-one-dimensional [(2+l)-d] physical systems, rather then the one-plus- 
one -dimensional [(l+l)-d] ones. For example, we will exploit the equivalence of NMM and 
the system of (2+l)-d Coulomb particles to compute correlation functions of the NMM in 
the scaling limit. We will also show that the partition function of NMM is a r-function of 
an integrable hierarchy containing (2+l)-d Burgers equations. 

The paper is organized as follows. To make the present exposition self-contained we 
reproduce in Section 2 the results of |J, devoted to the definition of the model, the derivation 
of the eigenvalue formula for its partition function and the integrability of the model. To 
emphasize the (2+l)-d nature of the NMM we notice following || that the partition function 
of the model can be interpreted as a classical partition function of Coulomb particles in a 
plane. In contrast, the partition function of the Hermitian matrix model (HMM) coincides 
with the classical partition function of Coulomb particles constrained to a string in a plane. 
We show that the NMM partition function can be written in a determinant form. This 
enables us to conclude that the partition function of the model is a r-function of Toda lattice 
with respect to holomorphic-antiholomorphic variations of the potential. 

In Section 3 we study the NMM with axially-symmetric potential. First we derive a 
determinant representation for the correlation functions of the model. We notice that the 
orthogonal polynomials associated with the integration measure are just powers of the com- 
plex variable z. This simplifies the analysis of the correlation functions considerably. We 
compute the correlation functions for the case when the matrix model potential is a mono- 
mial, V(M, M*) = (M ■ M^) k . The answer is expressed through degenerate hypergeometric 
functions. Thus NMM provides us with an example of a matrix model which is exactly (and 
explicitly) solvable beyond the Gaussian case even before taking the continuum limit. 

The two-point function of Gaussian NMM decreases exponentially at infinity, which should 
be compared to the power decay of the two-point function of HMM. While the latter can 
be obtained from quantum- mechanical computation in the system of free (1 + l)-d identical 
fermions in the external potential (see |JU| for details), the former follows from the analogous 
computation in the system of free (2 + l)-d identical fermions placed in the external magnetic 
field. 

In general we find that all correlation functions can be expressed in terms of a holomorphic 
function in one variable. This observation permits us to close a BBGKY chain of equations 
for the correlation functions and obtain a closed integro-differential equation for the two-point 
function. The existence of a single equation which determines the two-point function and 
thus, through determinant representation, all correlation functions, makes the NMM with 
axially symmetric potential exactly solvable. 

In Section 4 we solve the obtained equation in the continuum iV — > oo limit. The answer 
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for the two-point function is universal, i.e. potential-independent up to a change of scale. 

In Section 5, we develop a (2 + l)-d free fermion formalism and construct the free fermion 
representation of the partition function of NMM. This fact has its analogy in the theory of 
HMM, in which the partition function admits (1 + l)-d free fermion representation. 

In Section 6, we use the formalism developed in Section 5 to construct a family of com- 
pletely integrable hierarchies (labeled by an integer N) of non-linear differential equations. 
We show that the N = 1 hierarchy contains the (2 + l)-d Burgers equations (multidimen- 
sional generalizations of the original Burgers equation |§], were considered firstin f27| , see |l7j 
for a review). We thus call it the Burgers hierarchy. We give explicit solutions to it. These 
solutions generalize the Hopf-Cole solutions to the (2+l)-d Burgers equations, JTS] and 0. 

The family constructed provides a multidimensional extension of KP hierarchy which can 
be explained as follows. It is well known that the KP hierarchy can be formulated in terms 
of the pseudodifferential operator W(t, d) = 1 + t0i(t)9 _1 + w 2 (t)d~ 2 H ( [^2| and f43fl , see 



10 1 for a review). It is also known that KP hierarchy admits a reduction specified by the 
condition that W(t,d)d N is a dif ferential operator, see p4| , [05[| We call it the KP(N) 



hierarchy. The set of equations of KP hierarchy coincides with that of KP(N) hierarchy 
with TV — > oo. We show that the iV-th representative of our family of hierarchies of non- 
linear differential equations is a multidimensional extension of KP(N) hierarchy. Thus we 
call it the extended- KP(N) hierarchy. In the simplest N = 1 case the extended- KP(1), or 
equivalently, the Burgers hierarchy contains the (2 + l)-d Burgers equations, while the KP(1) 
hierarchy contains the (1 + l)-d Burgers equation. 

We then classify all formal solutions to the extended- KP(N) hierarchy. We find the one- 
to-one correspondence between the set of formal solutions and an open subset of an infinite- 
dimensional Grassmann manifold Gr (oo, N). This subset consists of all iV-dimensional sub- 
spaces of an infinite-dimensional complex linear space having a non-degenerate projection 
onto a fixed N- dimensional subspace. 

Next we show that the partition function of the N x N NMM is a r-function of the 
extended- KP(N) hierarchy. This is achieved by using the bosonization formula which can 
regarded as a (2 + l)-d generalization of (1 + l)-d bosonization formulae. 

Finally in Section 7 we discuss the Ward identities for the NMM. It is known that matrix 
model solutions to the KP hierarchy (like HMM) satisfy Virasoro constraints, i.e. they are 
annihilated by an infinite set of differential operators spanning a subalgebra of Virasoro 
algebra. This subalgebra is isomorphic to an algebra of infinitesimal holomorphic polynomial 
diffeomorphisms of the complex plane. We show that the partition function of NMM is 
a special solution to the extended- KP(N) hierarchy annihilated by the set of differential 
operators generating an algebra of all infinitesimal polynomial diffeomorphisms of the plane. 
The results of this Section can be used in the further analysis of the continuum limit of the 
NMM. 

The proofs of lemmas and theorems presented in the paper are placed between "<0" signs. 
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2 The Eigenvalue and the Determinant Forms of the 
Normal Matrix Model 



In this Section we will give the definition of the normal matrix model and derive the eigen- 
value and determinant formulae for the partition function. We will be utilizing the standard 
methods of the theory of matrix models (see for review). In the end of this Section, we 
will derive a connection between NMM and the Toda lattice hierarchy. 
The partition function of NMM has the following form: 

Z N =[ dti(T)e~ trV ^ M '\ (1) 

where F denotes the set of iV x N normal matrices, d[i(F) is a measure on F induced by flat 
metric on the space of all N x N complex matrices and V(z, z) is a function on C such that 
(HD exists. 

As usual, integral (|I|) can be reduced to the integral over eigenvalues {zi}f =l and {zi}f =l 
of matrices M and M< respectively. A corresponding calculation can be easily performed 
given the explicit expression for dfi(F) in the appropriate local coordinates on F. So let us 
sketch the calculation of d[i(F). 

First we notice, that any normal matrix M can be presented in the form M = UDU^ , 
where the D is diagonal matrix, U G U (N) is a diagonalizing matrix. This decomposition is 
unique up to multiplying U by a diagonal unitary matrix from the right. The corresponding 
equivalence class of U together with D defines a convenient coordinate system on F. Using 
these coordinates one can calculate the Riemannian metric induced on F by flat metric on 
the space of all matrices: 

|| 8M \\ 2 = trSM5M^ = tr(SD ■ 5D ] ) + 2tr(5u ■ D ■ Su ■ D ] - 5u ■ 5u ■ D ■ (2) 

where 5u = W6U is an invariant (Haar) length element on U(N). Due to the t/(A^)-invariance 
5u is well-defined in terms of our coordinates on F. Expressing || 5M \\ 2 through eigenvalues 
and matrix elements (Su)ij = (U'SU)ij we get 

N N 

|| 5M \\ 2 = ^SziSzi - SuijSujilzi - Zj\ 2 . (3) 

i=l *J=1 

On the other hand, || M || 2 = G ab l a l b , where {/ a } is a cumulative notation for the local 
coordinates on F, G a b is an induced metric on F. Then (jl(T) = det(GY^ 2 F[ a dl a (see e.g. |TT|| ). 
Combining these two formulae with (§) we obtain 

JV 

dfx{F) = dUlldz.dzilAiz)] 2 , (4) 

8=1 

where dU = Tli^j duijduji is Haar measure on U(N) and A (2) = det\z\~ ]i<i,j<N = 11^(^1 — 
Zj) is the Van der Monde determinant. 
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Finally, substituting (||) into flU) and integrating over the unitary group we arrive at the 
eigenvalue expression for Z N : 

N 

Zn = c(N) (n<M*e- V( * iA) )|A(z)| 2 , (5) 
JcN i=i 

where c(N) is the volume of the unitary group, a constant factor independent from V. From 
now on we replace it by unity. 

Next we will deduce a determinant formula for the integral (|5|) and relate the result to 
the solutions to Toda lattice hierarchy of differential equations. 

The reformulation of (||) in determinant form is standard and is based on the following 
formula: 

det[M lk ] ■ det[M jk ] = £ det[M iaU) M Mj) ], (6) 

aes N 

where M is any N x N matrix, a is an element of symmetric group S N . Applying (^) to the 
product of Van der Monde determinants A(z)A(z) in (|5|) we arrive at the desired form of the 
partition function NMM: 

Z N = N\ ■ det[Zij]i< it j< N , where (7) 
Zu = I dzdze- v ^z l - x z>- x . (8) 



Let us consider the potential V(z, z) in the form 

V t (z, z) = U(z, z) - J2(t k z k + t k z k ) (9) 



fc>0 



with tk = = tk for k 1. We are interested in the behavior of Zjv[14] with respect to 
variations of fs and t's. Note that Zjv[Vt] can be considered as a generating function of 
correlators in NMM with the potential U(z, z): 

<{trM\r -{trM^r---> NMM 
1 r N / 

— J CN Uidzidz^^^lAiz^i^trM^Y 1 ■ {trM^Y ' 



Z N Jc» f = \ 

~JN\yt\ \t,t=o 



Vi . (_\n... Z 



Zn dti 1 dti 2 

From (H) with potential given by (Qj we can easily deduce that 

dZ{j dZij , . 

-Qj- = z (i+k),j, -j^- = Z it{j+k) , i, j > 1. (10) 

Equation (0) together with ( [TOD mean that Z N [V t ] is an A^ th r-function of Toda lattice (with 
respect to complex variables t's and fs, see and |^7| for details). Note that the real version 



of the integral (|) together with potential (0) was studied in |33[], where it was referred to as 
a "scalar product model" . 

In Section 6 we will discuss a relation between NMM with an arbitrary potential and 
integrable systems which includes the one described above as a particular case. 
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3 Correlation Functions of the Normal Matrix Model 

In this Section we will analyze the structure of correlation functions of NMM with arbitrary 
axially symmetric polynomial potentials. We will see in this case that the two-point correla- 
tion function can be expressed in terms of the (analytical continuation of) one-point function, 
thus leading to the exact solvability of the model. 

It follows from (|J) that the distribution function of the eigenvalues in NMM with an 
axially symmetric potential is 

Pn(z,z) = ^\A(z)\ 2 e(-^ v ^), (11) 

where is the partition function of NMM. The n-point correlation functions (often called 
reduced distribution functions in statistical physics) are defined as follows: 

N\ r N 

RP (zi, ■ ■ ■ , z n ) = — / J] dzidziP N (z,z), (12) 

where the combinatorial prefactor accounts for the symmetry of the integrand. Note that 

R N ( z ) = (K z ~ zi)S(z - Zi)) N MM (13) 
and coincides therefore with the level density of the NMM. 



Correlation functions (12) can be presented in the determinant form. To prove it we note 



that the expression (|TT|) for the distribution function can be rewritten as 

P N {z, z) = —det[K(zi, Zj)]i<ij<N, ( 14 ) 

where 

N-l 

K N (zi,Zj) = 0k(z~i)M z j)i ( 15 ) 

A:=0 

and 

<t> m {z) = ^z m e-^ 2 \ (16) 

where c m > 0. The functions TO (z)'s introduced above are the orthogonal functions of the 
problem normalized by the following condition: 



dzdz4> m (z)(j) m i(z) = 8 m>m >, (17) 



which yields 



c 



7TC m / dx 



00 r 

x m e~ v ^ x > 



o 



1. (18) 
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Apparently, the orthogonal polynomials associated to the NMM with axially symmetric po- 
tential are just monomials in z. This explains the certain simplicity of the model. 

It is easy to check that the matrix K^(zi, Zj) is hermitian and it can be considered as the 
kernel of projection operator, i.e. the relation 



c 



dvdvK^{u,v)Kj<i{v,w) = Kn(u,w) (19) 
is satisfied. Therefore one can make use of the well-known result from the theory of matrix 



models (see | 30fl , page 89, theorem 5.2.1) to obtain the following representation of the n-point 



functions defined in (|T2j): 

RP(zi, ■ ■ ■ , z n ) = det[K N (zi, Zj)]i< id < n . (20) 
It is convenient to rewrite the expression for K N (zi,Zj) by substituting (|16|) into (|15D: 



2V-1 



K N ( Zi , Zj ) = £ c™ • & ■ ^eK-nNWfel 2 )) 

M^)e^ {N)2 - y(N)2 ). (21) 



Therefore, up to a factor which explicitly depends on the potential, Kn(z,w) is completely 
determined by the function of one complex variable k^{z ■ w). 

To illustrate the formulae obtained above let us consider the NMM with the Gaussian 
potential, V"( 1 2r| 2 ) = \z\ 2 , in the limit iV — > oo. A simple calculation gives the following 
answer for the kernel: 

K(z, w) = le^-^lMM 2 ) _ (22) 



7T 



Substituting (^) into ( p0|) we conclude that in the limit iV — > oo the one-point correlation 
function is constant and equals to -, whereas the 2-point correlation function in the same 
limit is 

R( 2 \z,w) = -e ( - |z " w|2) ]. (23) 

TV 2 

It follows from the last two results that the connected part of the 2-point function is equal to 
^2exp[—\z — w\ 2 } and decays exponentially at infinity. Let us note that the Gaussian NMM 
is equivalent to the Gaussian Complex matrix model (see P0f , Chapter 15). However this is 
no longer true for a non-Gaussian matrix model potential. 

Computability of NMM goes beyond the Gaussian case. For example it is possible to ob- 
tain a closed expression for K(z, w) when the NMM potential is a monomial in \z\ 2 : V(z, z) = 
2k , where k is a positive integer. Corresponding integral in ( |TBD is easily calculated to give 



\z 



the following expression for c m : 

k 



Trr(^i) 



m+TY' (24) 
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Substituting (H) into (|T]) one obtains the following; answer for the kernel at N = oo: 



^ (zw) p 



( p+1 

p=0 1 I J v 



k 



M e (-§N 2fc -|IH 2fc ) 



(25) 



where i^i,i(a; 6; z) is the degenerate hypergeometric function, the solution to the Kummer 
equation which is analytic at z = 0: 



d 2 uj .du 
z^^r + [p — z)— au = 0. 



dz 2 



dz 



(26) 



For k — 1, fl25|) reduces to fl2~2"|). It is well known (see e. g. [16[]) that the F lt i(a,b, z) has a 
Taylor expansion in z with an infinite radius of convergence and defines therefore an entire 
function on the complex plane. We use this remark below to analyze the correlation functions 
of NMM in the case of an arbitrary polynomial potential V. 

Suppose V(x) = Ylt=i tn,x n with td > 0. Then there are positive constants a\ and ci2 such 
that 



V(x) < ai + a 2 x d , x > 0. 



(27) 



This estimate can be used to prove that for any complex number u and a positive number R 
such that \u\ < R, 



k 



■N 



U 



for V(x)=Y, d n=1 tn*" 



< e 



^Ar(-R) for V(x)=a 2 x d 



(28) 



But the r.h.s. of (|2q) is nothing but a partial sum of a series converging by the remark 
above for any < R < oo. Therefore k(u) = limiv^oo k^ exists for any u such that \u\ < oo 
and defines an entire function on the complex plane. But an entire function is uniquely 
determined by its values on say positive part of the real line, which provides us with the 
connection between 2- and 1-point correlation functions. To make this connection explicit 
we notice that by definition 



^ (1) (M 2 ) 



k 



z? ] e 



-v(\z\ 2 ) 



(29) 



Denote by R^ (u) an entire function of a complex variable u which coincides with the 1-point 
function above for u = \z\ 2 . Such function exists as V is a polynomial and k(u) has been 
proved to be an entire function. Moreover it is unique, therefore is determined completely by 
the 1-point function. We finally notice that the connected part of the 2-point function can 
be expressed in turn in terms of R^(u): 



Rf\z,w) 



zw 



-[V(\z\ 2 )+V(\w\ 2 )-V(zw)-V(zw)] 



(30) 



which states the claimed property of the NMM. It's worth mentioning that there is the 
counterpart of (p0|) for any finite iV as well. 

Relation ( |3"0"D leads to an exact solvability of NMM. Let us discuss this point in some 
more details. It follows from (E0|) and (ETT) that all correlation functions of the model can 
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be expressed in terms of R^'(u). On the other side they can be considered as a solution 
to the certain chain of integro-differential equations which is discussed below. These two 
remarks permit one to derive a closed equation for the function R^>(u) which is what we 
mean claiming the exact solvability of NMM. To obtain such an equation we rewrite an 
expression for the distribution function ( JlT| ) m the following way: 

P N (z,z) = — eV ^* =1 11 11 ; ^>o 11 (31) 
Zn 



fl3T|) can be identified with the distribution function of classical two-dimensional Coulomb 
gas in equilibrium. The analogous equivalence has been widely exploited in the study of 
conventional matrix ensembles (see |3(J for review). 



Such an identification allows one to use the Liouville's Theorem of classical mechanics to 
derive a BBGKY chain of equations for the reduced distribution functions (|X2|) (see e.g. |32|j ). 
In particular we have the following equations, connecting one- and two-point functions: 

+ i?«(|,| 2 )^ = / dwdw R%\z,w)l z ln\z - w\, (32) 

and another one obtained from (|32|) by means of complex conjugation. Substituting ( p0|) into 
(j32|) we get a closed equation for R^(u): 



dR »}l Z?) + Ri\\ Z \ 2 )^ = i#(|*| 3 ) J dwdw R%{\w\*)^ z ln\z - w\ 

2 ' ' 



— / dwdw 



RV(zw) e-^ 1 ^^ 1 - 12 )-^-^]^^^-^! . (33) 

dz 



Equation (|3^), together with a reality condition R^^u) = R^(u) and normalization con- 
dition / dzdzRjy (\z\ 2 ) = N, can be used to determine the entire function R^\u) and thus 
all correlation functions of NMM. In what follows we will exploit ([33|) to study the iV —>■ oo 
scaling limit of NMM. 



4 NMM in the Continuum Limit 

The present Section is devoted to the NMM in iV — > oo scaling (or continuum) limit. We 
will find the corresponding limits of level density and the two-point correlation function. We 
will see that in the regions where the limiting level density is non-zero and smooth (no points 
of phase transition) the form of the two-point function is universal, i.e. independent of the 
NMM potential up to a rescaling . The scale itself is determined by the level density which 
happens to be inversely proportional to the square of correlation length. 

It is appropriate to mention that our search for universal answers in the NMM has been 
motivated by the study of universality in the hermitian matrix models, see e. g. ||. 

Consider the N x N normal matrix model with the potential N ■ V(\z\ 2 ), where V(x) is 
a real polynomial with positive coefficient in front of the monomial of the highest degree. 
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We assume V to be "convex," by which we mean that the following inequality holds for any 
z, w G C: 

V(\z\ 2 ) + V(\w\ 2 ) - V(zw) - V{zw) > 0, (34) 
and the equality in fl34|) is reached iff z = w. Setting w = z + e and expanding the l.h.s. of 



(34|) around e = we get the infinitesimal version of the convexity condition: 

AV(\z\ 2 )>0, (35) 

where A = dd is a two-dimensional Laplace operator. It is easy to check that any monomial 
potential of positive degree satisfies (0), therefore any polynomial of positive degree with 
positive coefficients does. 

We will be interested in the asymptotic properties of the model as N — > oo. We will 
consider two cases. First, we consider the limit N — > oo when the range of variation of \z\ 2 is 
of order 1, which corresponds to the limit of strong potential and large separation of levels. 
Second, we will investigate the case of strong potential and small separation of levels, i. e. 
the limit N — > oo when the range of variation of \z\ 2 is of order ~. 

For the first case we will seek the solution to ( j33|) in the form of the asymptotic expansion 
in the inverse powers of N, 



>(i) 

K aj — iv 

n=0 

The limiting one-point function is defined to be 



R^(u) = Nj:R^(u)N-\ (36) 



R^\u) = Rt\u)= lim (37) 

Substituting ([$6|) into (^) and equating the terms of the first order in iV we get the following 
equation for the limiting one-point function: 

R^(\z\ 2 )( 9V ^ 2) - j dwdwRS 1 \\w\ 2 )^ z ln\z-w\) = 0, 
which implies that 

(A - f dwdwR^(\w\ 2 )-^-ln\z-wl (38) 



dz J dz 
if we assume that i?W(|z| 2 ) ^ Q. To derive (|38D from (^) we have used the convexity 



condition and relation (pop. Relation (|30| ) was used to show that in the limit considered the 
connected part of the two-point function is zero almost everywhere. Note that we derived 
( |33"D from the Liouville theorem, so its direct consequence, equation (^), provides an exact 
description of the continuum limit of NMM. 

It is interesting to note that equation (|31f) can be interpreted as an extremum condition 



for the potential energy U[p] of the gas of charged particles of density p(|^| 2 ) = R^(\z\ 2 ) 
placed in the external potential: 

U[p] = J dzdzpV(\z\ 2 ) — - J J ' dzdzdwdwp{\z\ 2 )p(\w\ 2 )ln\z — w\ 2 . (39) 
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Such an interpretation was a starting point in the derivation of the level density in the 
hermitian matrix model presented in P0 |, however it required an assumption that the two- 



point function factorizes in iV — > oo limit. In contrast, in the NMM we can demonstrate the 
factorization of the two-point function almost everywhere by analyzing (PD|). 

Equation (|38| ) can be easily solved by differentiating both sides with respect to z. Noticing 
that ln\z — w\ is a Green function of two-dimensional Laplacian, i.e. Aln\z—w\ = nS^^z — w), 
we get 

i? (1) (M 2 ) = -AV(\z\ 2 ). (40) 
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So we conclude that the limiting level density (one-point function) of NMM at the point z 
is equal to ^AV^d^l 2 ), given that the limiting level density is not zero. This answer has a 
transparent physical meaning: the gas of charged particles which describes the continuous 
limit of NMM tends to screen completely the external potential. 

Having described the spectrum locally one can derive a complete picture assuming that 
the potential energy of the "gas" of NMM eigenvalues is minimal. In our case V is convex 
and the minimum energy assumption yields the following answer for the level density: 

R (1) (\z\ 2 ) = -AK(|^| 2 ) if |^| 2 < r (41) 

7T 

R {1 \\z\ 2 ) = 0if|z| 2 >r, (42) 
where the end point r of the spectrum is determined from the normalization condition 

dzdzR {l \\z\ 2 ) = 1. (43) 

\z\ 2 <r 

Consider now the case of small distances between eigenvalues, for which the two-point 
correlations are essential. Let x be a real positive number. We assume that the i? (1) (x) > 
and there is a closed neighborhood of x in the complex plane where the convergence of the 
r.h.s. of ( |3"TD is uniform (in physical terms, x is not a point of a phase transition in the 
N — > oo limit). We will study the limiting behavior of the two-point function of NMM 
in the vicinity of such a point. Note that the existence of points of uniform convergence 
does not contradict the information about the structure of correlation functions we have 
collected so far. For example, the uniform convergence of the sequence of entire functions 
R$ in the closed neighborhood U of x implies that the limit is a holomorphic function 
in this neighborhood. But according to (f40|) R^(x) admits an analytic continuation in 
some neighborhood of x given that R^(x) ^ (recall that V(x) is a polynomial function). 
Moreover, a direct verification shows that in the case of monomial potentials x is a point 
of uniform convergence if x is not too large and R^(x) ^ 0. We conclude therefore that 
the domain of our consideration is not empty. However, it does not consist of the whole 
plane: it is apparent that the endpoints of the spectrum defined by Q4"5| ) are excluded from 
our considerations, as the sequence of continuous functions cannot converge uniformly to a 
discontinuous one. Thus the results concerning the two-point functions presented below are 
valid only in the "bulk" of the spectrum. 
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Consider now the following holomorphic change of coordinates: 



z N = y/x- e c/ ^, w N = y/x- e v/v ^ \ u N = z N - w N . (44) 

The corresponding limit of the one-point function is R^(x, £ + rj) = lirriN^ 00 RN It is a 
matter of simple considerations to show that if liniAr^oo gN{u) = g(u) uniformly, u e U, then 
limAr^oo Qn{u ■ /at) = g(u) for any sequence {/n} converging to 1 and such that u ■ E U 

for any N. Using this remark we see that R^(x, £ + rj) = BP~'(x). Let us also remark that 

the expression R^'(x, C + rj) = R^'(x) solves the limit of the equation ( f44|) corresponding to 
small separations between levels. One can consider this remark as yet another consistency 

check of the assumption about the existence of the point of uniform convergence. Now we 

can substitute (44 into (|30|) and take the limit N — ► oo of fl3"0|) to obtain the following answer 
for the scaling limit of the two-point function: 

R^\x,C, V ) = [i?«(x)]V^ (1) WI^I 2 , (45) 

where we used the result ([401) and assumed that R^'(x) ^ 0. We conclude from ([45]) that 
the scaling limit of the two-point function of NMM exhibits universal behavior since the 
dependence on potential (which enters the answer through R^ 1 ' (x) only) can be eliminated by 
changing the scale and normalization. Notice that the square of correlation length is given by 
inverse level density. Equation (|45|) can also be interpreted as a completely universal relation 
between one- and two-point functions of the NMM. A similar phenomenon was previously 
observed in || for the hermitian matrix model and is referred to as "universality of the second 
type" . 



5 Free Fermion Representation of NMM 

The aim of the present Section is to rewrite the partition function ([!]) of NMM in the form 
of the free fermion correlator. This will permit us to describe an arbitrary variation of NMM 
potential by means of integrable system of differential equations and generalize therefore the 
results of Section 2. The technique we use is a straightforward generalization of free fermion 
methods developed in |4(| and reviewed in part in and their application to the theory of 



matrix models due to the ITEP group (see [33] for review). 



Consider an abelian group Z x Z consisting of pairs of integers. The corresponding group 
multiplication will be denoted by "+". To simplify notations we refer to this group as G. 
An abelian group Z acts on G as follows: 

G x Z -> G, 

— (i,j),m) i — > g • m — (i ■ m,j • m). (46) 



Alternatively, G can be presented in the following form: 

G = © meZ G m = G„ © G+, (47) 
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where 



G_ — © m <oG m ; G + — © m >oG m , (4£ 

and 



G w = g g G 



(l + m, -I) ; leZ , me Z. (49) 



Note that Go can be visualized as an antidiagonal of G = Z x Z. It is also worth mentioning 
that if g G G m and h G G n , then g + h G G n+m . Therefore, G is Z-graded and fl47|) is a 
decomposition of G into a sum of components of a fixed degree. 

Let us consider an infinite-dimensional Clifford algebra -4(G) over complex numbers which 
corresponds to G. We define it by means of the following set of generators and relations: 



-4(G) = ^1, i/j s , V> h ; g, h g G {if; s , ij h } = 8 S}h , |> g , ^ h } = = {^ g , ^ h } )■ (50) 

Let Tr be the right Fock space, an irreducible left -4(G)-module obtained by applying -4(G) 
to the right vacuum vector \vac), which is defined as follows: 

ij; s \vac) = 0, g G G_ ; ^{vac) = 0, h G G + . (51) 

In the same fashion we introduce the left Fock space JF L which is generated by the right 
action of -4(G) on the left vacuum vector (vac\ which is defined below: 

(vac\4> g = 0, g G G_ ; (vac\tp h = 0, h G G + . (52) 

There is a natural non-degenerate pairing between left and right Fock spaces, 

(vac\a, b\v ac)^j i— > (vac\a ■ b\vac), 

where a and b are elements of -4(G). This pairing is normalized as follows: 

(vac\l ■ l\vac) = 1. (53) 

From now on we will think of -4(G) as an algebra of linear transformations of or Tl and 
refer to elements of -4(G) as operators. 

To each element a G -4(G) one can assign its normal reordering with respect to the 
vacuum \vac), an element of -4(G), denoted as : a :. It is obtained from an element a by 
permuting the generators in each term of the sum constituting the element a in such a way 
that the fermion generators annihilating the right vacuum stand to the right from the rest of 
generators. 

A variety of vectors in the Fock space is provided by means of the following construction. 
To each finite set 5" C G we assign the following vector in Tr. 

\S) = JJ o s \vac), (54) 
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where o g = ip s if g G G + f] S and o s = -0 g if g G G_ f| Note that without fixing a linear 
order in the set Z x Z, the element (|53|) is defined up to a sign only. The same construction 
can be applied to obtain a vector in Tl which we will denote as (S\, 

(S\ = (vac\ n o g , . (55) 

where o g = ^ g if g G G + f] S and o g = ijj g if g G G_ f] S. It is easy to see that any element 
of Tr or Tl can be presented as a linear combination of vectors (|54]) or (|55D, respectively 

Consider a semigroup Q C G consisting of all pairs g = (m,n), where m and n are 
non-negative integers. We introduce the following element of the algebra -4(G): 

H{t) = E kJ s = E *g( E V^h +g ), (56) 
geQ geQ heG 

where 

tg = if cte#(g) » 0, (57) 

and we always suppose that to — 0. We will call H{t) a Hamiltonian operator or simply 
Hamiltonian. It is easy to check that the operators J g 's introduced in (|56|) commute, i.e. 
[J g , Jh] = J g ■ Jh — Jh ■ J s = 0, for g, h G Q \ {0}. It is also a matter of simple computation 
to verify that [J e ,ijJh\ = *ph-g and [J Sl tJjh\ = V'h+g- 

Before we introduced the Hamiltonian (|5"Bj ) all our considerations had not been different 
from the standard fermion construction as -4(G) is isomorphic to a standard Clifford algebra 
the generators of which are labeled by elements of Z, the isomorphism is being established 
by ordering the set Z x Z. The new development starts with the introduction of Hamiltonian 
( |56D since the presented commutation relations between J g 's and fermion generators depend 
on the group structure in G = Z x Z, but G is not isomorphic to Z as a group. 

It follows from the definition (^) that the Hamiltonian operator annihilates vacuum, 

H(t)\vac) = 0. (58) 

Define a ^-evolution of an element a G -4(G) as the following element of (the formal comple- 
tion of) A(G): 

a{t) = e H ^ae~ H ^. (59) 

Due to (|58| ) the t-evolution preserves the normal ordering, : a(t) :=: a : (t). It is not difficult 
to derive the following expressions for the t-evolution of fermionic generators of -4(G): 

Mt) = E PhWg-h, (eo) 

heQ 

Mt) = E Ph(-t)^ g+ h, (61) 

heQ 

where {p g } g eQ are generalized Schur polynomials, 

P g (t) = S Si0 + t g + — E + 4 51 t h rt h >rt h rn H . (62) 

h',h"eQ;h'+h"=g cJ - h',h",h"'eQ;h'+h"+h"'=g 
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The set {p g } g eQ defined above indeed consists of polynomials: the sum in the r .h. s. of (|6~2"D 
is finite as there is a finite number of ways to present an element of Q as a sum of elements 
of Q of positive degree. It is also worth mentioning that one can grade the polynomial 
ring C[th,h G Q] by setting deg(t h ) = deg(h). Under such assignment p s (t) becomes a 
homogeneous polynomial of degree equal to deg(g). 

We will also need a generating function for the Schur polynomials. Let z, z be complex 
coordinates in the plane. We set z g = z a ■ z b , where g = (a, b). Let also 

V(t, z)= E^ g - (63) 
geQ 



Then 



e 



V(t, z) 



£p g (^ g . (64) 



geQ 



The generating function ( pB| ) can be used to rewrite relations fl6"0| ) and ( KTj ) in a compact 
form. Let 

#1 = E^¥ g (65) 
gee 

be a free fermion field operator. Then a direct computation shows that 

e H ®4(z)e- H ® = e v <t*)i/>(Z). (66) 

To derive a free fermion representation of the partition function of NMM we will need a 
pair of operators called fermionic projectors: 

P + = : e E ^ G - ^ : , (67) 
P_ = : e -£««s + ^« : . ( 68 ) 

One can verify the following properties of fermionic projectors: 

P + ^ g = ^ g P + = 0, ifg€G_; (69) 
P-V>g = 4>sP- = 0, ifgeG + ; (70) 
P\ = P+ ; P 2 = P_. (71) 

Eq. ( [71] ) is a direct consequence of flBTSD and (0), which in turn result from the following 
equivalent representation of fermionic projectors: 

p + = n (i-^g), 

gec_ 

p = n (i-« g ). 

gec+ 

Finally we are ready to introduce the main object of our interest. Consider the following 
fermionic correlator: 

r(U,A(<S>)\t) = (U\e Hit) A($)\U), (72) 
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where U C Q is a finite subset of Q C G, 



1M*)= E = E ( 74 ) 

g eG+ g eG+ 

and $ is a real function (or a distribution) on C 2 . Our aim is to compute the correlation 
function (|72|) . Consider the following matrix of an infinite size: 

i($) gh = / dzdzdwdw §(z,w)w g f h , h, g G G + , 
Jc 2 

A($) g , h = «y g ,h, h, g G G (75) 

and A s ^(<&) = in all other cases. One can verify the following set of identities involving 
the matrix A($) g h: 

E A($) g , h ^ h A($) = A($W e , g G G. (76) 

heG 

To prove ([76]) we consider separately cases g G G_ and g G G + . In the former case 
relation (|76|) reduces to ^ g A($) = A($)?/> g , g G G_, which is true, as the operator A($) 
depends on fermionic generators labeled by elements of G + only The latter case follows in 
a straightforward way from the representation of the operator given below: 

00 i r 

A(Q>) — } — - / dz\dz\dw\dw\ ■ ■ ■ dz m dz m dw m dw m Q(zi, wi) ■ ■ ■ $(z m , w m ) 
"'■ J 

^+{z x ) • • -Mz m )P^ + {w^) ■ ■■^ + (w{~ 1 '- 1) ). (77) 
In what follows we will always require the non- degeneracy of the matrix A St b(<&) . We also 



wish to remark that operator A(<&) from ( [731) is a counterpart of the operator used in [22] to 
obtain a free fermion representation of conventional one-matrix models. 

Using relations ([76j) and the fact that A(&)\vac) = \vac) and applying Wick's theorem 
one can express r(U, A($)\t) in terms of two-point correlators: 

T(U,A($)\t)=det[(vac$ g (-t) E ^($Wkkac)] g)he£7 , (78) 

keG+ 

where the determinant is taken with respect to some linear ordering of U. Note however that 
the change of order is a unitary transformation which does not change the determinant. 

It is not difficult to compute two-point correlation functions entering (|78|). Substituting 
([75]) into ( [75|) and using relation ( ffify we arrive at the following expression for t(U, A(<&)\t): 

r(U,A($)\t) = det[Z g)h (t)] gih6C/ , (79) 
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where 

Z g , h (t) = J & dzdzdwdw z s w h $(f, w)e v{t ' s) (80) 

Suppose now that the set U C G and the function $ which parametrize the r-function are 
chosen to be the following: 



U 



N 



{gGG|g = (n,0),n = 0,l,---iV-l|, (81) 
$(z,w) = ® NM m{z,w) = 5(z - w) ■ 8{z - w) ■ e" |z|2 . (82) 
Substituting (|8l|) and fl32|) into (|79|) and (0), we see that 



Z itj (t)= f dzdz z l z> e£™.">° ( * (m '">~ 5m ' l5 "- l)2mr \ (83) 



which coincides (up to nonessential factor) with the determinant formula (f7|) and (||]) for the 
partition function of NMM with a polynomial potential equal to —V(t,z) + \z\ 2 . Thus we 
proved that 

Zn = j^t(u n , A($ NM m) 



t . (84) 



To conclude the discussion of the present Section we would like to mention that the 
partition functions of conventional matrix models (such as Hermitian, Unitary, Complex 
one-matrix models and Hermitian two-matrix models ) can be presented in the form ( |72|) 
under appropriate choices of U and $. To obtain the partition function of the hermitian 
matrix model, for example, set 



Then 



with 



$(z,w) = $hmm(z,w) = K z ~ w ) ■ 6(z-w) ■ 5{z-z) ■ e |z|2 . (85) 



t(U n ,A(<5> hmm ) 



t = cfe^jjcKijXAT-i (86) 



dx x i+i e *<»> oT ™ xm and T m = ]T t m) - 5 m>1 . 



Eq. (jS6|) is the determinant form of the partition function of hermitian one-matrix model 
with the potential — J2 P >o T p x p . 
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6 The Partition Function of NMM as a r-Function of 
the Extended- KP(N) Hierarchy 



In this Section we are going to derive a system of differential equations associated with 

In virtue of 



correlation function ([72). In virtue of (|83|) all results of the present Section apply as well to 
to the partition function of NMM with an arbitrary polynomial potential. 

The part of our construction dealing with free fermions relies heavily on methods de- 
veloped in § 



and 



sec 
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for a review. Our consequent analysis of the emerging 
hierarchy of differential equations shows that the original approach of and [53 to the 
theory of KP equations can be extended to the multidimensional well. 

Let Un and A be a subset of Q and an element of the Clifford algebra given by (|8ll) and 
(|73|) correspondingly (to simplify notations, from now on we denote ^4($) and A($) g .h as 
A and A gi h respectively). The following functions depending on z and fs are called wave 
functions: 



w p {z, t) = ± r S (87) 



(y N 




e H ®ip{z)A 


u N \ 




(u N 


e H{t) A 


u N x f 





(u N 






< 


u N 


e H(t) A 


U N 


) 



w q (z, t) 



where i/)(z) = SgeG^g-? g * s a f ree ne ^ operator; p = (1,0) G G, q = (0,1) G G. In our 
notations for the wave functions we suppress the dependence on iV and A which are supposed 
to be fixed. We assume that the common denominator in (|87|) and (|38D is not equal to zero 
when all t's are equal to 0. Then the wave functions make sense as formal power series in t's. 

There is a linear relation imposed on the wave functions which follows from the identity 
(|76|) . Let us explore it. Consider the Fourier decomposition of wave functions with respect 
to z: 



w p (z,t) = £ , w&t) = £ wl(t)z* (89) 

geG g eG 

It is easy to prove the following set of relations between coefficients of such decomposition: 

EAi,h<(t) = 0, £^g,h<(<)=0ifgGG_U^. (90) 
heG heG 

Here A s ^ is a matrix defined by (|75[). To prove (^) we observe that Wp(t) is proportional to 

(U N \iJN. P e H{t) ^ h A\U N }. 

Multiplying it by A gjh , summing over h, and using ( |76D we get (UN\i> P e H ® 'Aip g \UN), which is 
zero if g G G_ \JUn- The similar arguments apply if we replace w p (t) with w^(t). Therefore, 
( pl| ) is proven. 
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Relations ( p0[) have their counterpart in the theory of KP hierarchy (see e.g. |T(| ) and 
are of prime importance for our further considerations. Before one can make use of them 
however, it is desirable to rewrite equations (|90|) in the form of linear relations between the 
finite number of unknowns. This is possible because of the following representation of wave 
functions (E7i) and fl88|): 



ftt) = (z N -v+ E a s {t)z*)e v ^\ (91) 

g6f/jv 



Wci (z,t) = (z*+ E b s (t)z s )e v ^\ (92) 

where coefficients {a s (t),b s (t)} s£ u N depend on fs and the choice of A. It is easy to express 
them in terms of free fermion correlators, but we will not need the explicit expressions. To 
verify fl9T|) one can perform the following computation: commute the field operator ip(z) in 
the numerator of ( [37]) with e H( ^ using ( |66f) . Then notice that 

(U N \^ N . P ^ S = 0, g G G + \ (U N \J{N • p}) and ^ h e H ®A\U N ) = 0, h G G_, 



in which the first equality follows from the definition (|55|) of the generating vectors of Tl 
while the other one is a consequence of our choice of the operator A defined in (73) and the 
law (p0| ) of evolution of fermionic generators: 

ifj h e H ^A\U N ) = e H ^ E A^ gPg (-t)\U N ) = 0, if h G G . 

geQ 

It only remains to show that the coefficient in front of z N ' p in (|9~T|) is indeed 1. But according 
to ( p7|) this coefficient is equal to ^^§^w^j\fj^^ > which is 1. Representation (^Tj) has 
been derived. The derivation of (|9^) can be performed along the same lines. 

Using the fact that e v ^'^ = SgeGPgW-? s , where we assumed that p s = if g G G \ Q, 
and comparing representations of wave functions (^) and (fl9"T]), (|92"D ) one can rewrite the 
relations ([5t]) in terms of coefficients {a g , b g } sG u N - 

E \Ph- N . P (t) + E ak(*)Ph-k(t) I = if g G G (J U N , (93) 
heG \ kel/jv / 

E 4s,h (Ph-q(t) + E b k (t)Ph-u(t) I = if g G G |J U N . (94) 
heG \ k e^jv / 



Looking closer at relations (£)3|) and ( p4|) we see that the only non-trivial ones are those 
corresponding to g G Un- Indeed, if g G G_ then A g ^ = for h G G + in accordance 
with definition (|75|) . Thus left hand sides of relations (|93| ) and (|9~4l) are just identical zeros if 
geG_. 

Let us consider now the following pair of differential operators: 

W p (t,d) = d N *+ E a s (t)d s , (95) 

g6f/jv 

W^t,B) = d q + E ^W^ S . (96) 
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which we call wave operators. Here 3 g = (^-) 9P (^~) 9q an d g = <7 P ■ P + g q • q is an 
arbitrary element of G decomposed in terms of p and q. In what follows we will denote the 
complex variables t p and £ q as x and y, respectively. Thus W p and Wq are polynomials in 
^ and with coefficients depending on "times" t's. Evidently, 

«j p (t, z) = W p (t, 3)e v ^\ w n (t, 2) = W q (t, d)e v ^. (97) 



Noting also that d s ph(t) = Ph- S (t) for g G Q, we can rewrite relations ( fffi ) and ( p4|) in terms 
of wave operators, which will provide us with an alternative form of equations (|90|) suitable 
for our needs: 

W p h s (t) = = S G ^, (98) 

where h s (t) = J2h&G ^g,hPh(t)- Relations (^) can be viewed as two N x N systems of linear 
equations with respect to coefficients of wave operators. Here is an explicit solution: 

det(d°- p h, d l p h, ■ ■ • , d^'^h, d Np h, a (n+1) ' p K, ■ ■ • , d^-^h) , , 

a n-p(i) = =5 ~ ~ =3 , (99) 

P det{dP- p h,---,d^- x > p h) 

_ det0°- p h, d lp h, d^-Vvh, d^h, d^ n+1 > p h, d^-^h) 
n ' P ~ det(d°-Ph, ■ • ■ , d^-V-Ph) ' 

whre n = 0, 1, • • • N — 1 and h is a column vector with elements h g , g G Un- It follows from 
the results of the previous Section that the common denominator of ( p9|) and (|100|) is equal 
exactly to (Uj^\e H ^ A\Un) , and is therefore an invertible power series in £'s. Thus the wave 
operators are uniquely determined by conditions (j98D and their coefficients are formal power 
series in "times" . 

Let B be a ring of differential operators in ^ and with coefficients being formal series 
in variables t's. The statement bellow is a natural consequence of the formalism developed: 

Lemma 1 Let O G B be a differential operator such that 

Oh g = 0, g G U N . (101) 
Then there are differential operators b p G B and 6 q 6 B such that 

O = b p W p + 6 q H/ q , (102) 
in other words O = 0\ mo d w p ,Wq- Moreover one can choose b p to be of zeroth order in J^. 

The proof of Lemma 1 is given in the Appendix. The system of non-linear equations 
satisfied by coefficients of wave operators is a direct consequence of the formulated statement. 
To see this we will differentiate each of the relations ( |98|) with respect to £h and use the fact 
that -^~h s (t) = d h h s (t), a property of Schur polynomials. As a result we obtain the following 
set of identities: 

dWp ~Mr^ + Wp (q)(t ' ^ K = > if § G ^. (io3) 
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It also follows from (|98[) that 



[W p ,W q ]/i g = o, g eU N . 



dW p (q) (*,g) 
<9t h 



(104) 



+ W p (q) (t,d)0 h G B and 



Relations ( |103| ) and ( |104| ) state that operators 

[Wp,W q ] G B annihilate the set of functions h s , g G U^. Thus it follows from Lemma 
1 that 



dW p (t,d) 
dt h 



+ W p (t,d)d h = 



ot h 



, h G Q, 

i 

, h G Q, 



(105) 
(106) 



[W p ,W q \ = Y-W p = 



mod Wo 



(107) 



where F 6 B is completely determined by W p and W q . An explicit expression for the 
operator Y in terms of wave operators will be derived later and the answer is given in ( |129| ). 
Note that the r.h.s. of the ( 11071 ) is proportional to W p only, which follows from the fact that 
[W p , Wq\ is an operator in only 

For a fixed N, relations (|105|) and ( |106j ) constitute the system of non-linear differential 
equations for the unknown functions {a n (t), b n (t)}^Q subject to the constraint ( |107| ). More 
explicit form of this system will be presented later. 

At this point we must comment on the correctness of the definition ( |105| ) and ( |106| ) of 



equations of our hierarchy. Depending on how we mode out the parts proportional to W p 
and W q we can obtain seemingly different answers for the remainder. The reason for such 

Oh g = 0, g G Un^I consisting 

P ,W q , 



ambiguity is purely algebraic: the left B-ideal I = |0 G B 

of all differential operators annihilating h s with g G Un is not freely generated by W l 
there is the relation ( |107[ ) between generators. Moreover, this relation is the only one. In 
other words, we have 



Lemma 2 The left ideal I C 
relations: 



I 



B can be described as follows in terms of generators and 



[Wp, W q ] = 0\ mod Wf 



:io8) 



which means that for any O G I there are b p G B and 6 q G B such that O = b p W p + b q W Ci . 
Moreover, expressions & p Wp + 6qWq and ^Wp + ft^Wq are presentations of the same element 
of I if and only if there is an element c G B such that 

b p = b'-c- (W q + Y) and b n = b'+c- W p . 



Here Y G B is an operator defined in ( JUj ). A proof of Lemma 2 is given in Appendix. 
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Therefore we conclude that all possible ways to write the remainders in ( |105|) and ( |106| ) lead 
to the same answer if the relation ( |107| ) is taken into account; thus the the hierarchy of 
equations we care about is completely determined by (|105|) , (|106| ) and ( |107|) . 

The practical way of deducing differential equations from identities ( |1(J5 ), (106) and ( |1U7| ) 
can be extracted from the proof of Lemma 1 in the Appendix. To illustrate the result we 
present explicitly the simplest equations among ( |105| ), (|106|) and ( |1071 ) in the case when the 
set Un consists of one point, i.e. N — 1. The expressions for the wave operators in this case 
are 



The corresponding equations for h 



(2, 0), (1, 1) and (0, 2) can be written as follows: 
+ 2(d x a)a = d^a, 
+ 2{d x b)a = d 2 x b; 



db 



dt 



(2,0) 



da 

d£(i,i) 
db 

<9£(1,1) 



+ (d x a)b + (d y a)a = d x d y a, 

+ (d X b)b+(dyb)b 



d x d y b; 



da 

<9£(0,2) 

db 

(0,2) 



dt 



+ 2{d y a)b = d y a, 

2(dyb)b = d'b. 



The additional condition (|107|) takes the form 

d x b = 0. 



d y a 



(109) 
(110) 

(111) 
(112) 
(113) 
(114) 
(115) 
(116) 

(117) 



Let us rewrite equations (111) - (|116|) and ( 117 ) in terms of real variables and real- valued 
unknown functions. It follows from the reality condition imposed on the potential V(t, z) 
that t( 



'-(2,0) 



£(0,2), £(i,i) 



t, 



',(1,1) and x = £(1,0) = £(o,i) = V- Thus a(t) = b(t) and we introduce 
the following new real variables: 

Vl {t) = a{t) + b{t) , iv 2 (t) = b(t) - a(t), 
r 1 = x + y , ir 2 = x — y 
n = £(i,i) , £(2,o) = 2(r 2 + iT 3 ),£(2,o) = 2(r 2 - ir 3 ). 



Note that as far as transformation properties are concerned, Vi(t) with i — 1, 2, is a covector 
field (one- form). Relation ( |117| ) states that this one- form is closed: 



d (vi(t)dr 



0. 



:il8) 
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Equations (|111|) - ( 116| ) written in terms of Vi(r, r) acquire the form 



+ (v a ■ V)(4) = A g J a , a = 1, 2, 3, (119) 



where V = is a gradient operator, g a ^s are metric tensors, (gi)ij = 5 i j 1 {g2)ij = 

i(cr 3 )jj, (gs)ij = (<Ji)ij and 03 = diag(l, —1), ai = antidiag(l, 1). Matrices cr's are the Pauli 
matrices. An operator A 9a = g^didj is a Laplace operator of two-dimensional space equipped 

with metric g a . Finally, v a {i) = ( v^t), v^(t) ) is a vector field corresponding to a covector 



field yvi(t),v 2 (t) j in the presence of the metric g a , v l a = g^Vj. Note that g x is a Euclidean 

metric, while gi and g% are equivalent Minkowski metrics. 

Relations ( |119|) for a = 1 (the Euclidean case) are called two-dimensional Burgers equa- 
tions, f27j. As a result we see that two-dimensional Burgers equations are included in an 



infinite hierarchy of non-linear differential equations, equations ( |105| )- flT06| ) and ( 1 07| ) for 



N = I. This hierarchy is completely integrable in the sense that we can find all solutions in 
the class of formal power series, see Theorem 1 below.| 

Burgers equations in one, two and three spatial dimensions together with continuity 
equation and potentiality condition ( |118|) can be used to model potential turbulence. We refer 



the reader to [117) for a review and original references. This book also describes an application 
of three-dimensional analogue of equations ( |119|) to the study of large scale structure of the 
Universe. The integrable structure we have discussed can prove useful in the systematic 
analysis of the development of the turbulence in the models based on Burgers equation: 
the integrable structure of the Burgers hierarchy implies that the dynamics of the system is 
constrained to the invariant subspaces of the phase space (or " state space" in the terminology 
of J2BJ). The addition of the small perturbation in the form of the random force (see e.g. 



46| ) destroys integrability and the system moves towards chaos through the deterioration 
of invariant subspaces in accordance with Kolmogorov-Arnold-Moser theory. The presented 
picture is very close to the existing scenarios of the development of the turbulence ( Hopf- 
Landau scenario for example, see [2£| for details) but has a chance to admit a complete 



qualitative treatment. Moreover, the integrability (in the sense of the presence of integrals of 



motion) plays an important role in the description of fully developed turbulence of [38] and 



For example, Polyakov's anomaly introduced in is an anomaly of the conservation 
law. We hope to investigate the consequences of complete integrability of Burgers hierarchy 
for the theory of Burgers turbulence in the near future. 



In the meantime let us discuss the solutions to ( |119|) subject to the condition ( |118| ). 



Solutions to ( |118] ) which are defined at every point of the (r 1 , r 2 )-plane are of the following 
form: 

v i (t)=d j *(t), (120) 

where \I/(t) is an arbitrary function of t's. It is called a potential of the covector field Vi(t). 
Comparing (|120|) with (|99]) and ( |100| ) we see that the whole class of solutions to equations 

2 By the way, the original (i.e. one-dimensional) Burgers equation [|J also can be included in the integrable 
hierarchy constituting a certain reduction of KP hierarchy This statement can be easily extracted from the 
results of |35|. 
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119 ) is given by 

= ln(U x \e H ® A\U X ) = In^J dzdz[J dwdw$(z, w)]e v ^|. (121) 



We see that ty(t) is a generalization of the Hopf-Cole (|]T3, [0]) solution to the Burgers 



equations ( |119| ). The quantity inside the square brackets in the r.h.s. of the above equation 
is determined by initial conditions. The choice of $ from ( |73|) in accordance with fl82|) makes 
it clear that the (1 x 1) NMM solves the (2+l)-dimensional Burgers hierarchy. We also see 
that the Hopf-Cole substitution 

Vi(t) = d^(t) = d l ln{r{t)) (122) 

linearizes the (2+l)-d Burgers hierarchy, (|105|) and ( |106| ) at N = 1, subject to the constraint 
( |118| ). Note that the function r appearing in ( |122j ) can be considered as a generating function 



for the solutions of the (2+l)-d Burgers hierarchy and plays in this sense a role of the so 
called r-function of this hierarchy, see [fEJ where the notion of the r-function was introduced. 



Now we wish to analyze in some details the structure of equations ( |105| ), (|106|) and (|107|) 
for an arbitrary N. Their equivalent form is the following: 



dW p (t,d) 



+ W p (t, = 0* p W p + Oj^Wq, h e Q, (123) 



at 



h 



9W lf 8) + W^t, d)^ = 0\ p W p + OlJV^ he Q, (124) 



[W p ,W q \=YW p , (125) 

where Y and Oij with i, j = p, q are elements of B and according to Lemma 1 we can choose 
these operators in such a way that Oj iP with i = p and q are of zeroth order in J^. Such a 
choice permits us to express the r.h.s. of ( |123| )-( |125|) in terms of wave operators alone. To 



obtain such an expression we have to define first the right inverses of the wave operators. A 
way to do this is the following. We set 



W p 1 = d~ N Y,dnd- n , 

n>0 



n>0 



where {c? n },^ are formal power series in t's, while {e n } c £ =0 are differential operators in d x 
with coefficients being formal power series in t's. Operators {e n }^l and series {<i n }^L are 
uniquely determined by equations 

W p -W- 1 = l = Wq-W~ 1 . (126) 

The multiplication operation "•" in the above equations is defined by the Leibnitz rule. 
Relations ( |126p considered as equations with respect to unknown quantities {d n }'^ =0 and 
{e n }^L possess a unique solution. 
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Multiplying equations (|123 ) , and (|124|) by W 1 from the right, extracting the differential 



parts and using the fact that the O h 's are already differential operators and the Of p s, i = p 
and q, are of zeroth order with respect to d y we find 

Oj q = (wi^W" 1 ) , i = p and q, (127) 

where the subscript "plus" denotes the operation of extracting the differential part of an 
operator, i.e. if 0(t, d) = S g eG c g9 g then (o(t, d)^j = Z) g eQ c gd g - 

To compute operators O pp and 0^ p we substitute ( |127| ) back into ( |123j ) and ( |124j ). 
Multiplying the resulting equations with W~ x and projecting them onto differential operators 
in d x only we get: 

Ot = ((W pM d h W^ 1 )_W Ci W- 1 ) ,i = p, q (128) 

v / (+,o) 

where the subscripts denote the following operations: for any pseudodifferential operator O, 
O- = O — + and 0(+,o) is a projection of O onto differential operators in d x . 

It remains to compute the operator Y entering the r.h.s. of ( |125|) . This differential 
operator is of zeroth order in d y and is therefore equal to 

Y=([W p ,W (l ]W; 1 ) {+fi) . (129) 



Substituting ( |127| ) - ( [L29|) into ( |123|) - ( |125|) we obtain a system of equations governing the 
evolution of the wave operators: 

+ Wp & = ((W p frW- x )-W«W- x ) {+fi) W p + (Wpfrw-^+W^ g e Q (130) 

+ w ^ = ({W^W-^-W^W; 1 )^ W p + (if/iy, 1 )^,, g e Q (131) 



[W P ,W n ] = ([W^W- 1 ) W p . (132) 



Equations ( P-30|) and (|131|) together with the constraint (|132|) constitute a hierarchy of non- 
linear differential equations which can be viewed as a generalization of the Sato equations in 
the theory of KP equations (see for review). Burgers equations ( |119| ) with the condition 



(|118|) give the simplest examples of equations ( |130| ) - (|132|) for N = 1. 



For large enough N equations ( |130| ) and ( |131j ) contain first n equations of KP hierarchy, 
n « N . To see this, consider the equation (|130|) for W p when g = n ■ p with n > 1. In 
order to present the answer in the standard form we set W p = W and t n = t n . p and obtain 

dW 

+ Wd n x = (WdgW-^+W, n > 1. (133) 
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The set of equations (|133|) constitutes a certain reduction of KP hierarchy, which we describe 



as follows. Take the solution to KP hierarchy ( (J10|| ) , 

W = 1 + w 1 (t)d~ 1 + w 2 (t)d- 2 + ■ ■ ■ , (134) 

which satisfies an additional condition of W ■ being a differential operator. Then the 
operator W = W ■ solves equations ( |133| ). We call the hierarchy of equations ( |133| ) the 
KP(N) hierarchy. This hierarchy has been described in details in [|3l| . 



The KP hierarchy itself can be viewed as a limit of KP(N) as TV tends to infinity. Such 
a limit makes sense due to the stabilization of equations in KP(N) hierarchy: its n-th 
equation is independent from N if n << N. However the hierarchy ( |133| ) at finite N can 
be of independent interest as well. For instance, the exact version of the statement of the 
footnote 2 is that KP(1) is an integrable system containing the (l+l)-d Burgers equation. 
We can interpret equations ( |130|) and ( |131|) with (|132|) as an integrable extension of 



KP(N) hierarchy (|133|) to higher dimensions. Therefore, we call such system the extended- 
KP(N) hierarchy. The N = 1 example considered above supports such an interpretation: 
the extended- KP(1) hierarchy, or equivalently the (2+l)-d Burgers hierarchy, is a natural 
extension of the KP(1) hierarchy, or equivalently the (l+l)-d Burgers hierarchy. 

Our main result about the extended- -fT.P(iV) hierarchy is that it is a completely integrable 
extension of the KP(N) hierarchy, i.e. all solutions to the extended-i ; CP(A^) hierarchy are 
of form (p9|) and (|100|). To be precise we have the following theorem. 



Theorem 1 Let C[[t]] be a ring of formal power series in t s ,g 6 Q ; with complex coeffi- 
cients. Let B be a ring of differential operators in d x and d y , with coefficients belonging to 
C[[t]]. Operators W p and W q e B of the form 

N 

W p (t, d x , d y ) a n (t)d n x ~\ (135) 

71=1 

N 

W tl {t,d x ,d y )=d y + J2K{t)d:-\ (136) 

n=l 



where a n (t) and b n (t) 6 C[[t]] for n — 0, • • • , N — 1, solve the system ( \13Q ) - (\133j ) if and 
only if there exists a set {h n (t)}^ =1 of N elements o/C[[t]] such that 

The Wronskian W x (hx, ■ • ■ , hjsr)(i) of h\, - ■ ■ , /ijy 
with respect to the variable x is an invertible element o/C[[t]], (137) 

d s h n (t) = ^-h n , n=l,---,N; g G Q, (138) 
dt g 

and 

W p h n (t)=0 = W (l h n (t), n=l r --,N, (139) 
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Note that (i) operators W p and Wq defined in ( |135|) and ( |136| ) are the explicit form of 
the wave operators defined in and (^); (ii) the condition ( |137| ) implies in particular 
the linear independence of hi(t), ■ ■ ■ , /iAr(t) G C[[t]]; (iii) relations ( |139|) constitute a system 
of linear algebraic equations for the coefficients ai(t), ■ ■ ■ , a^(t) and bi(t), • • • , &jv(£)- This 
system has a unique solution due to the linear independence of functions hi(t), • • • , hjy(t). 
Thus Theorem 1 indeed describes all solutions to ( |130|) - ( |132| ) and the explicit form of these 
solutions is given by (^) and (|100| ). 

The proof of Theorem 1 is based on the following lemma. 



Lemma 3 Let W p and Wq be any two elements of B of the form ( 135 ) and (\13(\ ) satisfying 



the relation ( \133j ). Then dim(KerW p f]KerWq) = N over the ring C[[t p ,t q , t]] of formal 
power series depending on {t s | g G Q \ {p,q}}- Moreover the basis in KerW p f]KerWq 
can be chosen to satisfy condition &137{ ). 

The proof of Lemma 3 is presented in Appendix. However the proof of Theorem 1 based on 
Lemma 3 is so short and transparent that we present it here. 

^> Suppose that operators W p and Wq G B solve the hierarchy ( |130|) - ( |132j ). Let us 



fix a basis h\(t), ■ ■ ■ , h]y(t) of KerW p f] KerW^, which satisfies condition ( |137|) . Lemma 3 
states that any element of C[[t]] annihilated by W p and Wq can be decomposed into linear 
combination of hi(t), ■ ■ ■ , h N {t) with coefficients depending on all £ g 's, g G Q \ {p, q}. Let 
us apply operator equalities ( |130| ) and ( |131| ) to • • • , hjqit). The r.h.s.'s of the results 



are identically 0; substracting = a ^ f '^^ from the l.h.s.'s, we get: 

W p D s hi(t) = = WqD s k(t), g G Q,i = 1, ■ ■ ■ N, (140) 

where D s = ^ d s . Therefore D s hi(t) G KerW p C\KerWq for each i and g. Then there 

exist N x N matrices A s (t) with g G Q independent of x and y such that 

Djait) = XK(*)]#i(*)- (141) 

3 

These matrices are not unrelated. The fact that [D s , D^] = together with the linear 
independence of basic elements hi(t), ■ ■ ■ ,h>N(t) yields 

r)A f) A 

~di ~ + ^] = 0, g, h G Q \ {p, q}. (142) 

These zero-curvature-like conditions imply that there is an (x, y)-independent non-degenerate 
matrix B(t), such that 

A S = d ^B-K (!43) 

Non- degeneracy of B(t) means that detB(t) is an invertible element of C[[t]]. Consider a 
new basis of KerW p f] KerWq defined by 

^) = ^B|(t)^(t), i = l,---,N. (144) 
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Substituting (fT43|) and ([144]) into ( |i~41~D we see that 



D g hi(t) = 0, i — 1, ■ ■ ■ , N. (145) 



Therefore the condition (|138|) of the theorem is satisfied by elements h±(t), ■ ■ ■ , hi(t). The 



condition ( pL37| ) is satisfied as well, since W x (hi, • • • , hpf)(t) = detB{t) ■ W x (hi, • • • , /ijv)(t), 
and the r.h.s. of this relation is invertible in C[[t]]. Thus we have proven the "if part of 
the theorem (for each pair W p and Wq solving ( |130| ) - ( |132| ), there exists a set of linearly 



independent functions hi(t), • ■ ■ , /ijv(^) satisfying conditions (|137|) - ( |138| )). 

From our considerations which led to the hierarchy ( |130| ) - ( |132| ) we know that the inverse 
statement is also true: any pair of operators W p and Wq annihilating N functions that satisfy 
condition ( |137| ) and ( |138j ) solves the extended- KP(N) hierarchy. This concludes the proof 
of Theorem 1. <0 

Theorem 1 implies a geometric description of the space of solutions to the extended- 
KP(N) hierarchy. Before we can give such a description some additional notations are to be 
introduced. Let L C C[[t]} be a complex linear space consisting of elements of C[[t}} which are 
annihilated by operators D g = — d g , g G Q. Consider a set of all N- dimensional linear 
subspaces of L which we identify with an infinite-dimensional Grassmann manifold Gr(oo, N). 
Let us remind that GV(oo, N) is defined as a set of all iV-dimensional linear subspaces of C°° 



(understood as a Tychonoff product), see [|13|] for details. We equip Gr(oo,N) with the 
structure of topological space by declaring that the set £> 7 consisting of iV-dimensional linear 
subspaces of L having non-degenerate projections on the finite-dimensional linear subspace 
7 C C[[t]] is open. It follows from standard theorems of analysis (see e. g. [ffijfl ) that the 
set B = {i3 7 } 7C c[[t]j together with an constitutes a base for the topology on Gr(oo,N). 

Let ^GV(oo, iV)^ be an open subset of Gr(oo,N) consisting of all A^-dimensional linear 

subspaces of L having a non-degenerate projection on the subspace II of C[[t]) spanned by 
{x n }n~Q. The space of solutions to the extended- ATP (AQ hierarchy can be now described as 
follows. 



Corollary 1 There is a one-to-one correspondence between the set of solutions to ( |i 3(\) 
( [73^ ) and points of ( GV(oo, N) 



<) First of all let us prove that two sets of functions hi(t), ■ ■ ■ , /ijv(t) and hx(t), ■ ■ ■ , hisr(t) 
both satisfying ( |137| ) and ( |138[ ) determine the same pair of operators W p and Wq from 



and ( |136|) which annihilate them if and only if these two sets of functions are related by a 
constant non-degenerate linear transformation. 

The "if statement is a direct consequence of the Kramer's formula (see (p9|) and ( |100|) ). 
Conversely, suppose that W p and Wq annihilate two sets of functions hi(t), ■ ■ ■ , h^it) and 
hi(t), • • • , h^it) satisfying ( |137| ) and ( |138| ). Each of these sets span KerW p f| KerWq. There- 



fore, there is a non-degenerate N x N matrix M(t) independent of x, y such that hi(t) = 
J2j M- (t)hj(t). Applying the operator D s to both sides of the last equality and using the 
fact that D g hi(t) = D g hi(t) = 0,i = 1,---,N, and the linear independence of elements 
h±, • ■ ■ , fix, we get ^® = 0,g G Q. Therefore, M is a constant non-degenerate matrix. 
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Thus any two sets of N elements of C[[t]] obeying (|137|) and ( pL38| ) satisfy ( |139| ) iff they are 
related by a constant non-degenerate linear transformation. 

Theorem 1 implies that we have just constructed a one-to-one map from the space of 
solutions to the extended-XP(A^) hierarchy to the set of N- dimensional linear subspaces of 
L or, equivalently, GV(oo, N). This map is not onto: a linear subspace of L belongs to the 
image iff one can find a basis {h n }^ =1 of this subspace such that ( [1371 ) is satisfied. But 
this condition is equivalent to the non-degeneracy of the projection of our subspace onto the 
subspace II described above. Thus the image of the map in question is exactly (GV(oo, iV)) 
and the Corollary is proved. <0 

Finally, let us discuss a relation between the solutions (p9|) and ( |100| ) to the extended- 
KP(N) hierarchy and the partition function of NMM or, more generally, the fermionic cor- 
relators ([72|) with U = Ujy. Consider the so called vertex operator: 

X(t, z) = e v ^e~ V( "^' * ~ P) , (146) 

where V(~ ■ gf—,z~ p ) = Y^nyo^dT—z ~ n ' v . It is proved in the Appendix that the wave- 
function w p which encodes half of the solution to the extended-KP(iV) hierarchy can be 
expressed through the fermionic correlator: 

This equation can be viewed as a generalization of previously known one- dimensional bosoni- 
zation formulae JlO| . Knowing the fermionic correlator (|72|) we can compute w p (z,t) or, 



equivalently the set of functions {a n (t)}^ = Q. The other half of the solution is given by the 
set {&„(*) }n=o X - 

Even though we do not have at the moment explicit formulae for the functions 6(t)'s in 
terms of the fermionic correlators (^) we believe that functions 6(t)'s are also determined 
by t(Un, A,t). The reason for such a conjecture is very simple: all solutions at hand are 
determined completely by the matrix A g h with g G Un and h G Q which in turn can be 
restored from t(Un, A, t). This conjecture has been verified in case N — 1, see ( |120| ). So 
we conclude that fermionic correlators t(Un, A, t) with operators A from (|73| ) play a role of 
r-functions for the solutions (|9~9D to the extended- i'i'P(iV) hierarchy. 

We hope to continue the investigation of the structure of hierarchies ( |13U| ) - ( |132| ). One of 
the most important questions to be answered here is following: is there a universal integrable 
system from which one can obtain all hierarchies ( |130| ) and ( |131| ) for iV = 1, 2, • • • by means 
of reductions? The answer to this question is not clear. The reason is that equations (|130| ) 



and ( |131| ) of the extended- KP(N) hierarchy do not stabilize as N becomes larger. This is 
readily seen from the fact that the number of elements of Q of a given degree d grows with 
d. Thus the relation between extended-KP(A^) hierarchies and this hypothetical universal 
structure should be different from the known relation between KP(N) hierarchies and KP 
hierarchy. 
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7 Ward identities in the NMM 



It is well-known (see [31], [15|, [19], see [33] for review) that partition functions of hermitian, 



unitary, complex, etc., matrix models exhibit invariance with respect to a subalgebra of 
an algebra of holomorphic diffeomorphisms of a complex plane, or Virasoro algebra with 
zero central charge. This invariance can be presented in the form of the so called Virasoro 
constraints imposed on the partition function of a matrix model. It is also known that 

]2"T| , see |23j for a 



Virasoro constraints can be rewritten in the form of loop equations, 
review. These equations happen to be exactly solvable in certain scaling limits thus providing 
a powerful tool for a study of matrix models. 

The aim of the present Section is to demonstrate that the partition function of NMM 
is also subject to an infinite set of constraints. These constraints generate a subalgebra of 
Diff(C) algebra of all infinitesimal diffeomorphisms of the complex plane. 

We start with some heuristic considerations intended to unveil the reasons for the ap- 
pearance of a subalgebra of the Dif f(C) algebra in the normal matrix model. 

Consider a family of O-dimensional field theories with action V parametrized by the set 
of coupling constants {tki}k,l>o, where t kii = t^ k : 



N 



v = Y. E 

i=l kd>0 



(148) 



This family is equivalent in the field-theoretical framework to the NMM itself as for each 
fixed set of coupling constants ( 148|) gives a NMM potential entering (||). The algebra of the 
following reparametrizations acts in the space of field theories ( |148j ) : 



Z$ ~\~ 6 Za Za 
~ I p ~m+l n 



(149) 
(150) 



where m,n > 0. These reparametrizations can be presented as vector fields on the parameter 
space of NMM with the potential (|T45[): 



s m , n v 



e w m , n V + e w m , n V, 



where 



m,n 



ki 



dt hA 



Wr 



o 



ki 



k,l>0 ^ L k+m,l+n k[ > dt k+n j +m 

These vector fields obey the following commutation relations: 



(151) 



(152) 



(jp — m)w m+Pjn+q , 
(jp — m)"u; m+ p jn _|_ 9 , 



in which we recognize the relations for the Dif f(G) algebra generated by operators z 



(153) 



n+l ^ m _d_ 



Z 



dz 



and z n+ z n -^=. Operators ( |152| ) span a subalgebra of Diff(C) consisting of all polynomial 
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reparametrization of the plane. So we conclude that the subalgebra of Diff(C) describes 
classical symmetries of NMM. 3 

Having discovered the Dif f{G) symmetry on the "classical" level we gain the hope that it 
translates to the "quantum" level in the form of corresponding constraints (Ward identities) 
on the partition function ([![). To see this directly let us perform the change of variables given 
by (|149|) , (|150|) in the integral (|). As a consequence of the fact that the integral doesn't 



depend on the choice of the integration variables we obtain the following set of identities: 

jv ; / 

/ f[dz i dz i \A(z)\ 2 e- tr ^<>o tklMkM - J2 pt pq tr(M m+p M) n+q 

i=l \ p,q>0 

1 1 m 

+-(m + l)tr(M m M n ) + - £ tr(M m ~ p M n )trM p 

2 2 p=0 

+ £ E E c x ~ v n zr 1 -^ = o, (154) 

/ p=0 k^l Zk Zl J 

and its complex conjugate. Here M and M are diagonal matrices with Ma = Zi,Ma = z^ 
Our goal is to rewrite ( |154j ) in the form of differential constraints applied to the partition 



function itself. The obvious obstacle for doing so is the term with double summation in the 
left hand side of ( |154j ). To overcome it we present a fraction z z k Z z z \ i n the following form: 

* k ~ * l - {z k - Zl f r due-^- z ^ +r >\ (155) 



where rj > is an infinitesimal parameter. Note that if Zk = Zi the right hand side of ( |155| ) is 
equal to - = 0, therefore (|155|) can be considered as a continuation of the fraction at hands 
to all values of The replacement of ^Ef 5 with the integral ( |155| ) in (|154|) alters the 

integrand on the set of measure and doesn't change the integral itself. Using this remark 
we can substitute ( |155| ) into (|154|) to arrive at the desired result: 



W mn Z N = ; W mn Z N = 0, (156) 

where m, n > and 

d (I, _ d 



W m ,n = J2 kt kl^l + {-~(™+!) 



dt r 



k,l>0 ^ L m+k,n+l 

1 A d d r , ,,./ 8 ,.n d 



+ oE^T aT + oEt ^ e 



2 p=0 ^m—p,n dtp Q 2 JO \<9t(),n+l— p ^m+l,p 

+ -r^e-^ 2 » «_)}, ( 157 ) 

TF m , B = + {_I( m + i)JL_ 



3 The reasonings above are inspired by [pT| 
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1 ™ d d l ^ /-oo / d A «9 

+ ~X)a a + _ E / da;e X ( q e a 



Here 



D =5 5-55-55 + 5 J, (159) 
and 5 and 5 are the shift operators defined as follows: 

d d 

" at ' (160) 

- d d 

6{ d^^dt ' (161) 



The arrows above operators 5, 5 in ( |159| ) indicate the direction in which these operators apply. 
The operators e _tjD in ( |157| ) and ( |158| ) act within the brackets encompassing them. Note that 



the first terms in the expressions ( |157| ) and (|158| ) for the operators W mn and W mn coincide 



with the "classical" Diff(C) generators (|152|). The additional "anomalous" terms in the 



curly brackets of expressions for W mn and W mn appear from the variation of the measure in 
the integral (|5|). 

One can check directly that operators ( |157|) and ( |158| ) satisfy the commutation relations 



(|1 53|) , thus providing us with the Diff(C) constraints for the NMM. It follows from (|153|) 



that operators W m fl and W n fl generate two commuting copies of the subalgebra of Virasoro 
algebra. Thus the partition function of NMM model obeys Virasoro constraints as well. There 
is a nice feature of Virasoro generators W m $ and W n $ following from their representation 
(|157|) and (|158|) : the tedious terms involving integrals over uj are absent in this case, so the 



expressions for Virasoro generators W m; o and W n fl are similar to the ones appearing in the 
conventional matrix models (cf. p3[): 



r\ -| r\ -| 777, r\ r\ 

W m ,o = E k t kl - -(m + 1)— - + - E (162) 

k,l>0 ui m+k,l z ui m,0 z p= Q (Ji m-p,0 <JL pfl 

w n , = E i t»-J- - V n + + \ £ df—^- (163) 

k,l>0 ol k,n+l z L/C0,m z p= Q <J l 0,n-p ^0,p 

The generators of Diff(C) given by ( |149| ) and (|1 50|) can be presented in the form z n l m 
and Wrn,n = z n l m , where l n ,l m generate two commuting copies of Virasoro algebra. In 
this representation the invariance with respect to both Virasoro algebras implies Diff(C)- 
invariance. One might hope that, correspondingly, the invariance of the partition function 
with respect to W mi o and W nt0 of equations ( |162|) and (|163|) implies the invariance with 
respect to Diff(C) generated by ( |157| ) and ( |158j ). Presently we do not have any results 



which could justify such a hope. 

In conclusion of the present Section we would like to make the following comment. Usually 
the derivation of Ward identities in the conventional matrix models deals with the partition 
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function written as an integral over the subset of the set of all matrices. Unfortunately, such 
a representation of the partition function of the NMM is unknown due to the non-linearity 
of the space of normal matrices. However, we have managed to derive -Di//(C)-constraints 
for NMM starting from its the eigenvalue form (H). Actually, this approach has its own 
advantages. For example we are now able to apply similar considerations to analyze the 
structure of Ward identities in the models which are genuinely non-matrix. 



For example, consider a one-dimensional (3-model, which was originally introduced in [12 



see |$0| for a review. This and similar models reappeared in recent literature in quite different 



contexts such as study of chaotic systems [35] and fractional statistics [33 
function of the model is 



2~4| . The partition 



z N (p,t) 



N 

Y[dxi\A(x) 

i=i 



k=l 



(164) 



where (3 is any non-negative number and the integration goes over R^. The orthogonal, 
hermitian and quaternionic matrix models are the particular cases of the /3-model for the 
values of (3 equal to 1,2 and 4 correspondingly. In general ( |164j ) can be thought of as a 
classical partition of the two-dimensional Coulomb beads threaded on a string; (3 is identified 
in this case with the charge of a bead. The expression in the exponent has a meaning of 
the external potential. Exploiting the invariance of the integral ( |164|) under the change of 
variables xi — > Xj + ex™ +1 , n > —1, we obtain the following set of Virasoro constraints satisfied 
by the model: 



l£z N (/3,t)=0,n>-l, 



(165) 



2 t^ n dtk dt n _k 



fz = KJV k, yyv n—lt " ^"n k=l ^fc+ra 

To our best knowledge this simple observation was not done before. 



(166) 



8 Conclusion 

We see that NMM is a highly tractable yet nontrivial model possessing rich structures. The 
correlation functions of NMM with an axially symmetric potential can be expressed in terms 
of a holomorphic function of one variable. This holomorphic property leads to the universality 
of the correlation functions in the scaling limit. In the particular case of monomial potentials 
this holomorphic function is expressible in terms of degenerate hypergeometric functions. The 
NMM admits a free fermion representation. Using it we find that the behavior of the partition 
function of the NMM with respect to an arbitrary variation of the potential is governed by 
a completely integrable system of non-linear differential equations. This integrable system 
constitutes a multidimensional extension of the KP(N) hierarchy. In the simplest case when 
N = 1 it contains (2+l)-d Burgers equations. The partition function of the NMM is subject 
to .Di//(C)-constraints which reflect the symmetries of the model. 
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From the results of this paper we can foresee the folllowing possible generalizations and/or 
developments of. For example there is an integrable hierarchy corresponding to an appropri- 
ate subset U C Q which is different from the set Un used in this paper. By "appropriate" we 
mean that the ideal of the ring B consisting of the operators annihilating the set h g (t), g G U 
is finitely generated. An example of such a set is served by {g G Q,\deg(g) < n}. 

Another way to generalize our construction is to start with an arbitrary finitely generated 
Z-graded abelian group G and fix a semigroup Q C G in such a way that the Z-gradation 
of G induces Z + -gradation of Q. Let us take for example G = ZxZxZ, Q = {(m, n.p) G 
G | m,n,p > 0}. In the simplest case analogous to the case N = 1 considered in the 
paper we will have three wave operators Wi = di + a,i(t),i = 1,2,3 subject to constraints 
[Wj,Wj] = 0. These constraints have a solution cij(t) = d^(t) and one of the simplest 
equations of corresponding hierarchy of equations written in terms of the potential \I/ reads 

<)X] ' + {d x *){d v d z V) + (d y V)(d x d z V) + (d z V)(d x d y V) 



= (d x d y d z *) + (d x *)(d x V)(d x V). 



This equation has the following solution: 



* = In 



which can be verified by direct substitution. 

It is worth mentioning that integrable systems in 2 + 1 dimensions which are different 
from KP hierarchy have been investigated recently using Lax formalism (see e. g. [|I| and 
p8| ): generalized Schroedinger equations, Davey- Stewartson equation (a generalization of 
non-linear Schroedinger equation to (2+1) dimensions) and their close relatives. It will be 
interesting to study the Lax structure of 2 + 1-dimensional Burgers hierarchy ( extended- 
KP{1)) introduced in Section 6. 

The following question seems to be extremely interesting to explore: is there a possible 
connection between NMM and multidimensional gravity analogous to the well-known re- 
lations between HMM and two-dimensional gravity? To address this question we need to 
formulate an analog of Feynman rules for the computation of the coefficients of the asymp- 
totic expansion of the partition function of NMM in the case when coupling constants are 
large. Corresponding analysis is complicated by the non-linearity of the space of normal 
matrices, but there are no principle obstructions for the progress. 

The asymptotic expansion of the partition function of the hermitian matrix model (with a 
non-polynomial potential) carries topological information about the moduli space of Riemann 
surfaces |3?| . What can be said in this respect about the partition function of NMM? 

As far as algebraic geometry is concerned it will be interesting to study an analog of 
quasiperiodic solutions to KP hierarchy fZB\ in the context of extended- KP(N) hierarchies. 



Such solutions could correspond to algebraic surfaces on the one hand and to points of 
Gr(N, oo) on the other, thus providing us with a sort of multidimensional Krichever con- 
struction, 
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9 Appendix 

9.1 The proof of Lemma 1 

Recall that wave operator W p defined by (|95|) is of order N in d x while W q from (|96"D is of 
the first order in d y . Any operator O G B can be presented in the form 

Mi M 2 

^ = E X! O mil m 2 d'^' 1 d y n2 , 
mi=0 rri2 =0 

where coefficients C's are elements of C [[£]]. Therefore, 

Mi 

O = ^2 ( ^m 1 ,M2^T 1 ^ 2 + terms of lower order in d y . 

mi 

Or, 

Mi 

O = ^2 Cmi,M 2 ^r i ^ 2_1 W / q + terms of lower order in d y . 

mi 

So, applying induction we see that 

O = b q W q + O', 

where O' G B is of zeroth order in d y . If the degree of O' with respect to d x is less than iV 
we are done. If not, we have O' = Ck9 x + lower order terms = Cxd x x ~ N W p + lower order 
terms. If the order of the omitted terms is greater than iV we continue the process. After a 
finite number of steps we will arrive at the following representation of the operator O': 

Of = b p W p + J] c s (t)d «, 

where b p is a differential operators in d x only. Therefore an operator O can be presented as 
follows: 

o = b p w p + & q iy q + J2 c s (t)d g - ( 167 ) 

gel/jv 

Applying equality ( |167| ) to the function h h (t) with h G Un we see that 

E c g (t)(a g /i h ) = o, heu N . (168) 

But ( |16^ ) is a system of linear algebraic homogeneous equations with respect to {c s } ge u N 
with the determinant which is an invertible formal power series. Thus c g = 0,g G Un and 
Lemma 1 is proved. 
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9.2 The proof of Lemma 2 

First of all let us introduce the grading of the ring B by setting deg(d x ) = 1, deg(d y ) = N. 
The degree of an element of B is then by definition the maximal degree of monomials in the 
linear combination constituting this element. We will prove Lemma 2 basing on the following 
proposition: 

Lemma 4 Any element O G B admits a unique representation of the form 

0= E Cn^W^WZ 1 , (169) 

711 ,712 >0 

where deg(c nitn2 ) < N and c nii „ 2 = if n\ » or n 2 >> 0. 

We will prove Lemma 4 later. Now let us show how the statement of Lemma 2 can be deduced 
from Lemma 4. 

Consider the following sequence of maps: 

0^BAB©B^I^0, (170) 



where a(b) = ( -6(W q + Y),bW p ),b G B and (3(c u c 2 ) = Cl W p + c 2 Wq, (d,c 2 ) G B © B. 

Lemma 2 is equivalent actually to the statement of exactness of this sequence. So, let us 
check the exactness of (|170|) . 

First, we see that the map a is monomorphism. Really, a(b) = implies in particular 
that bW p = which yields b = 0, as the ring B is an integral domain. Thus, Ker(a) = 0. 
Second, the map (3 is epimorphism. Really, it follows from Lemma 1 that for any O G I there 
are ci, c 2 G B such that O = c\W p + c 2 Wq = 0(ci, c 2 ). Thus, (3 is onto. 

Finally, we have to verify that Im(a) = Ker(f3). It is easy to see that Im(a) C Ker((3). 

Really, ■ a(b) = p((-b ■ (W q + Y),b- W p j) = b ■ ([W p , W q ] - YW P ), which is zero in 



virtue of the relation (|132[) . Let us prove the opposite inclusion, Ker{f3) C Im(a). If 
(ci,c 2 ) G Ker((3) then for any b G B (c 1 ,c 2 ) - a(b) = (ex + 6(W q + Y),c 2 - bW p J G 
Ker((3). By Lemma 4, c 2 = E nii „ 2 >o Qn^W^W^. Choose 6 = E m ,n 2 >o Qn 1+ i,n 2 W^W^ . 
Then, (ci,c 2 ) - a(6) = (<n + 6(W q + ^),E„>o ?o,n^ q n ) G Ker{(3). Again, by Lemma 4, 

Ci + b{W^ + Y) = En u n 2 >0 Pn u n 2 W^W^. So, 

E Pm,n 2 W™W™ +1 + E ?o,»C = o 

ni,7i2>0 n>0 

But the l.h.s. of the equation above is a representation of in the form (|169|) , it follows then 
from the uniqueness of such representation that p ni ,n 2 — for ni,n 2 > and q^ n = for 
n > 0. Going back we conclude that (ci, c 2 ) — a(6) = 0, which proves that Ker(j3) G Im(a). 

It remains to prove Lemma 4 which states the existence and uniqueness of decomposition 
(|169 ). Take any element O G B of order less than (n+ 1) • N with n > 0. It can be presented 
in the following form: 

71 

o = E + (171) 

fc=0 
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where the degree of q k is less than N and the degree of O' is less than k ■ N. Thus, 

0-0' = q WZ + q (d; - W$ + ■ ■ ■ + q n d n x N 

= qoW^ + q[d n y - l d^ + ■■■ + q' n d n x N + D u (172) 

where deg(q' k ) < N and deg(Di) < n ■ N . 

Suppose that we have proved that for some m such that < m < n, 

0-0' = c WZ + ■■■ + c m WZ- m W™ + d m+1 d n y - m - l d x m+1 ^ N + ■■■ + d n d n x N + D m+1 , (173) 

where degrees of operators q and dj are less than TV and deg(D m+ i) < N. Then we see that 

O-O' = c W^ + ■■■ + c m W^ m W™- N + c m+1 ^ q n - m - 1 Wj m+1 > w 

+ rf m+1 (a;- m - 1 ^ +1 > Ar - w ^- m - l w^ m+1 > N ^ + ■■■ + d n d n x N + D m+1 

= c W^ + ■■■ + c m+1 W q "- m - 1 W^ m+1 > 7V + d' m+2 d^~ m - 2 d x m+2 > N 

n®x N + An+2, 



+ ■■■ + d'd n j N + D 



where degrees of operators q and d'j are less than N and deg(D m+ i) < N. At some point 
we defined c m+ i = d m+1 . Thus we proved by induction, the base of which is ( |172|) and the 



induction hypothesis is ( |173| ), that 

n 

= J2c k W*W£- k + 0" 1 (174) 



fc=0 



where deg{0") < n ■ N . Note that if deg(0) < N, then the existence of ( |169| ) is clear. So, 
using ( p.74j ) to generate induction in degree we verify the existence of representation (|169|) 
for any O G B. 

To prove the uniqueness of (|169|) we must show that the equality 

E w^X" ( 175 ) 

ni,n2>0 

implies c„ 1>n2 = for all ni,n 2 > 0. We assume of course that deg{c ni , n2 ) < N . Suppose that 
n ■ N < deg(J2 ni ,n 2 >o c ni ,n 2 W^W^) < (n + 1) ■ N with n > 0. Then' " 

W n 2W ^ = c nun2 W^ 2 W^ + O'" = 0, with deg(0"') <n-N. 

ni+n 2 =n 

Thus c 0j ndy + ( terms of lesser degree in d y ) = 0. Therefore, c 0in = 0. Continuing step-by-step 
considerations of terms of the highest degree we see that c nii „ 2 = for n\ + n 2 = n. But this 
contradicts the assumption that deg(J2m,n2>o c «.i,n 2 

WqW™ 1 ) > n ■ N with n > 0. Thus all 
coefficients c nii „ 2 in the l.h.s. of ( |175| ) for rii + n 2 > are equal to and we are left with the 
equality cq,o = 0. The uniqueness of ( |169|) is thus proved and so is Lemma 4. 
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9.3 The proof of Lemma 3 



First we will construct N linearly independent elements of C[[t]] annihilated by W p and W q 
and satisfying condition (|137| ). Then we will prove that any element of C[[t]] annihilated by 
W p and W n is a linear combination of these elements with coefficients independent of x, y. 

Let hi,---,hN be iV linearly independent elements of C[[t]] generating KerW p . They 
can be chosen to have the following form: hf. = x k ~ 1 + ( higher order terms in x) with k = 
1, • • • , N, so that the condition (|137|) is satisfied. Applying (|132|) to these elements, we see 



that Wqhi(t) G KerW p . Therefore, there exists an N x N matrix E(t) independent from x 
such that 

W q h % {t) = Y J E{t))h,{t) (176) 

3 

Consider an invertible matrix F(t) solving the equation 

VWI-^^W- (177) 



Such a solution always exists in the formal category: it is easy to verify that F(t) = I + 
F 1 {t)-y + F 2 {t)-y 2 + ■■■, where/ is an N x N identity matrix and F m (t), m > are N x N 
matricies independent of y. They are determined one-by-one from the recursion relation 
F m+1 (t) = ^ ri Y, t+j=m E l {t) ■ Fj(t) with m > 0, E{t) = Ei>o^W^ and F (t) = /. The 
evaluation of our solution on complex numbers for y can be written as usual in the form of 
path-ordered exponent. 

Consider now a set of functions /ii,---,/ijv defined from: hi(t) = ^2AF{t)]\hj{t), i = 
1, • • • , N. Then W^h^t) = 0, i = 1, • • • , N due to ( [1751) and the fact that the matrix F(t) 
solves equation ( |177| ); and W p hi(t) — 0, i = 1, • • • , iV due to the fact that F(t) is independent 



from x. Thus hi, • ■ ■ , /ijv £ KerW p f] KerW^. Due to the non-degeneracy of the matrix F(t) 
and the fact that the elements hi, ■ ■ ■ , /ijv were chosen to satisfy the condition ( |137| ) it follows 
that elements hi, ■ ■ • , /ijy also satisfy the condition (|137|) 



Moreover, these elements generate the intersection of the kernels: suppose there is f(t) e 
C[[t]] such that f(t) G KerW-p^KerW^. Then f(t) G KerW p as well. Therefore, there are 
coefficients di(t), ■ ■ -(^(t) independent of x such that 

f(t) = J2d l (t)h l (t). (178) 

i 

Thus all we have to prove is that coefficients d{S are independent from y. Applying Wq to 

( |1 78| ) we obtain that Y^,id y (di(t)jhi(t) = 0. Then by linear independence of h^s over formal 

power series depending on y and t g , g G Q \ p, q we conclude d y (di(t)) = 0,i = 1, - ■ ■ ,N. 
Lemma 3 is proved. 



9.4 The proof of bosonization relation 



Here we present the proof of ( |147| ) 
P_(-t) = e- H{ $P_e H ®. Following pO 



First, consider operators P + (—t) = e H wp + e H W and 
we will call them Clifford operators. The identities 
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below || express the properties of the Clifford operators which will be important for us: 

(vac\P+(-t) = (vac\e H(t) = (vac\P-(-t). (179) 

The proof of ( |1 79| ) is not completely straightforward and was not explained in ||. There- 
fore let us outline the proof. 
Consider the state 



(vac 



e H Wp + (-t). (180) 



We wish to compute £(e|. Note that H(t) = e- H ^H{t)e H ^ = E g eG+ t g J g (-et), where we 
used the definition (EBI) and the agreement (1591). Therefore 



^-(vac\e H ^H(t)P + (-t) 
de 

= E h E (vac\e H ^M-^ s+ h(-et)P + (-t) 

g eG+ heG+ 

= E h E (vac\^ s+h e H ^P + (-t) = 0, 

g eG+ heG+ 

where the first equality is due to the properties (|69|)-([71|) of projector operators and the 
last equality follows from (|52|) . Therefore the state (e| is in fact independent of e. Finally, 
observing that (vac\P + = (vac\ we obtain the desired result: 

(vac\e m = (vac\e m P+(-t) = (e = 1| = (e = 0| = (vac\P + (-t). 

The first equality in ( |179| ) is proved. The second one can be verified along the same lines. 
Next we need to prove that 

(U N \e H ^^(z)AP + = X(t, z)z ^ N - 1 > p (U N \^ l) . p e H ^AP + . (181) 

The proof is the following. Using ( |6^) we can perform the following computation: 

[P + {-t)U N (-t),iP{z)) AP + 

= e v ^e~ H ^[P + U N ,i;(z)}e H ^AP + 

= e y M e -B®2 (N-i)v £ r^[P + U N , ^ N -x- n y P ]e H{t) AP + 

n>0 

= e^e~^.^- 1 >P[P4-^)^(-^),V-i).p]e^ ) AP + 
= X(t, z)e~ H ^z (^- 1 )-P[P + (- /U )C/ JV (- )U ), ^ {N _ iyp }e H ^AP + , 

where Y\ S GU N i J s = Un and /x g = - if g = n ■ p with n > and /i g = otherwise. A 
calculation shows that 

P + (-^=P + -e( g ~ PE ^^ s+p{ -" ) ) 
Using this result and observing that iph.P+ = 0, h G G_ we find that 

[P + (-t)U N (-t),iP(z)]AP + = X(t,£)i'^ 1 )-P e -^P + [^(-/i), ij {N ^ iyp }e H ^AP + . (182) 
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Applying the operator equality ( |182|) to the left vacuum and using the property ( |179| ) of 
Clifford operators we see that 

(vac\U N e H(t) ij(z)AP + = X(t, z)^' 1 ^ ' (vac\U N i) {N „ iyv e H ® AP + , (183) 

where we used that (vac\Ujsr(— fj,)ip^N-i)- P = {vac\U N^{N-iy-p-, which can be verified directly. 
Taking into account that (vo,c\Un = (Un\, we obtain ( |181| ). 

Finally let us apply ( |181| ) to \Un) and use the fact that P + \U~n) = \Un). Dividing the 
result by r(UN,A,t) we arrive at the equation ( |147| ). 
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